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The present work proposes the concept of induced percolation over multiple-object systems, so 
that features such as the number of merged clusters can be used as a relevant measurement. The 
suggested approach involves the expansion of the objects while monitoring the evolving clusters. The 
potential of the proposed methodology for characterizing the spatial interaction and distribution 
between several objects is illustrated with respect to synthetic and real data. 
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Great part of the current interest in statistical physics 
and complex systems targets critical phenomena, such as 
percolation p). Characterized by abrupt changes of the 
system properties as some parameter is varied, critical 
phenomena are of interest because of their rich and com- 
plex dynamics. As a matter of fact, while the calculation 
of critical exponents for continuous phase transitions has 
been identified as one of the greatest achievements of 
theoretical physics over the last 25 years @, the study 
of geometric percolation in graphs has provided one of 
the main motivations of the new area of complex net- 
works 3] . Much interest has been focused on the analy- 
sis, modeling and simulation of natural critical phenom- 
ena, whereas relatively little attention has been given to 
artificially-induced critical phenomena, namely the en- 
forcement of the objects in the system to undergo some 
imposed dynamics — e.g. expansion or dilation — in 
terms of a control parameter, so that critical phenomena 
can be induced. Although the intrinsic relevance of such 
investigations is potentially limited by the imposed evo- 
lution dynamics, it is argued and illustrated in the cur- 
rent work that such experiments present great potential 
for the characterization of the spatial distribution of ob- 
jects while taking into account their shape and topology. 
This perspective provides a novel approach to the old 
problem of geometric characterization of multiple-object 
physical systems, which has been approached from sev- 
eral perspectives including densities (e.g. additive 
functional (e.g. Jp), fractal dimension (e.g. 0,13) an< ^ 
lacunarity (e.g. |8|). The importance of obtaining an ef- 
fective methodology for the characterization of the spa- 
tial distribution and interaction in such systems is greatly 
enhanced by the myriad of related physical and biolog- 
ical systems, including statistical physics, astrophysics, 
geology, material science and even genetics (e.g. the 
spatial patterns of gene expression Q). Although the 
proposed methodology can be immediately extended to 
continuous spaces, the following developments are lim- 
ited to discrete structures as these are naturally required 
by several computational implementations and applica- 
tions. In this article we start by presenting the proposed 
induced-percolation concepts and methods and proceed 



by validating and illustrating them with respect to uni- 
form point densities, perturbed hexagonal lattices and 
real data regarding the spatial distribution of retinal gan- 
glion cells (e.g. Q). 

Let fi be the A''— dimensional Euclidean space where 
the multiple-object system is embedded and S be the set 
containing the coordinates (x, y) of each point of each 
object, in any order. The definition of the subsets of 
points corresponding to each object can be done by con- 
sidering the respective point connectivity or by imposing 
previous knowledge about their classes (in this case a 
single object may include more than one connected com- 
ponent). In this work, every pair of points a and b such 
that \x a — Xb\ + \y a — yb\ = 1 is linked through an undi- 
rected edge, so that each object corresponds to a con- 
nected component (or cluster) in fi. Let also R £ f2 be 
a set delimitating the region of £1 where the artificially- 
induced dynamics is allowed. The Euclidean D-ring of 
the objects in O is obtained by determining the set of 
points of R which are no further away than D from the 
set S [l| • Only the distance values di allowed in the or- 
thogonal lattice, i.e. di = \J x 2 + y 2 for some x, y belong- 
ing to f2, are considered in order to obtain full precision 
in spatially discrete spaces ^3 • Observe that the density 
of such possible distances tends to increase with the dis- 
tance. The dilation of S by the distance di can now be 
defined as the set of points obtained by the union of the 
original object and the respective D = e^-ring. Figure ^ 
shows a multiple-object system (a), the dilation of such 
a system for distances up to D = 10 (b), which can also 
be understood in terms of the exact distance transform 
of the objects 01 • 

The percolation dynamics considered in the present 
work involves the dilation, up to a maximum distance 
c?m, of the objects by each subsequent exact distance di, 
a scheme henceforth called exact dilations |l0j | . The exact 
distance di, or its respective index i 0], can both be con- 
sidered as percolation parameter. Two are the immedi- 
ate consequences of such an enforced dilation procedure: 

(i) the area of the objects increases monotonically; and 

(ii) contacts between objects occur at specific distance 
values, inducing the respective merging of the involved 
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FIG. 1: The dilations of co-existing uniform particle densi- 
ties (a) and the respective number of mergings C(s) and h(s) 
function. 



clusters and possible critical phenomena. By associating 
each object to a node and incorporating edges between 
pairs of objects whenever they first touch one another, 
it is possible not only to merge the objects into larger 
clusters, but also to obtain the hierarchy induced by this 
agglomerative procedure. 

Several measurements, such as the size of the maxi- 
mum cluster 0, , could be considered in order to iden- 
tify possible critical phenomena while the parameter i 
is increased, but the current work focuses on the cumu- 
lative number of clusters merged at each successive dis- 
tances G?i, which is henceforth expressed as C(di). Wc 
consider a square window with area A. In order to pro- 
vide spatial-scale uniformity, we make the change of vari- 
able C(ki = log(di)). Now, considering uniform distribu- 
tion of points (i.e. Poisson) with density 7 as a reference 
model, the rate of change at kj, while going from C(kj) 
to C(fcj+i), can be maximized by the finite difference 



q(kj) — P(kj)/Akj, where P(kj) — -fA is the expected 
number of points inside R and Akj = fcj+i — kj. In order 
to have the same change q(ki) for any value of ki while 
considering uniform point densities inside the fixed win- 
dow, we choose the changing rate at any specific value 
of kj (in this article j = 1) as a reference and make 
the normalization as given in Equation ^ By adopt- 
ing the minimal nearest neighbor distance as the proto- 
type distance for the uniform distribution with density 
7, we can use the mean expected value of such a dis- 
tance, given as di = exp(ki) = l/(2^j) ^jjj m order 
to relate 7 and fej. By using this expression, we have 
that P(ki) = 7A = A/(4d 2 ) and P(kj) = A/(4dj 2 ), so 
that P(kj)/P(ki) — exp(2ki — 2kj) and it is now pos- 
sible to rewrite Equation ^ a s Equation Now, the 
derivative of the number of cluster mergings in terms 
of ki for any type of point distribution can be ap- 
proximated as in Equation |3| which is used henceforth. 
This equation already takes into account the fact that 
kj—i = 0. An interpolated version of w(fcj) can be ob- 
tained by first embedding it onto the continuous space s, 
such that s = ki at i, as expressed in Equation and 
subsequently using the Parzen interpolating scheme [T^ 
given in Equation [5] to obtain the function h(s), where 
g a (s) — 1/ ^/2tt/ <jexp{— 0.5(s/a) 2 ) is the normal density 
with standard deviation a and '*' stands for convolution. 
Figures ^b) and (c) illustrates the distributions C(s) 
and h(s), respectively, for the multiple-object system in 
(a). The two co-existing Poisson densities in Figurc^a), 
which were generated with 71 = 0.05 and 72 = 0.002, 
are clearly identified by the two peaks obtained in func- 
tion h(s). Each of these peaks have respective maximum 
value at si = 0.728 and S2 = 2.673, yielding estimated 
densities 71 =1/(4* exp{2s 1 )) = 0.058 and 72 = 0.0012. 
Considering the interferences between the two co-existing 
points distribution and the relatively poor sampling of 
the loss dense distribution by the adopted 400 x 400 points 
window, such results can be considered as being reason- 
ably accurate. Moreover, the fact that only a fraction of 
the denser distribution present in the window is duly ex- 
pressed by the lower peak obtained for that case. Observe 
that such a multimodal distribution could by no means 
be obtained by using nearest-neighbor distance statistics. 
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In order to further validate the proposed methodol- 
ogy, it has been applied to a series of uniform point dis- 
tributions with increasing densities 7 = l/(4eP), where 
d = 1, 2, . . . , 15 are the characteristic spatial scales, esti- 
mated by nearest-neighbor distances. A 200 x 200 points 
window was adopted. The cases obtained for d = 3, 7, 11 
and 15 are shown in Figure [3 Two interesting results 
have been identified: (i) a reasonably good estimation of 
the typical densities can be obtained by using the sug- 
gested methodology; and (ii) the functions h(s) practi- 
cally do not change their shape or heights for the dif- 
ferent cases. It is also interesting to observe that the 
obtained functions are similar, but not identical, to the 
nearest- neighbor distributions given by Equation El also 
presented in Figure [3 

The potential of the proposed methodology for quan- 
tifying the degree of disorder Q in point distributions 
was also evaluated by considering hexagonal lattices with 
increasing levels of uniformly distributed noise. More 
specifically, the basic hexagonal distribution was charac- 
terized by basic cells with sides of 10 points, with pro- 
gressive perturbations added to the points coordinates. 
A 200 x 200 window was also considered for this case. 
It is clear from the obtained results that the perturba- 
tion degree is immediately reflected by the dispersion of 
the obtained functions h(s), while the typical point den- 
sity (indicated by the position of the peaks) is kept fixed. 
The sharpest peak is obtained for the noiseless hexagonal 
structure. The nearest-neighbor distributions for respec- 
tive density is also shown (solid line) superposed to the 
h(s), whose relatively narrow dispersion indicates that 
the point system is far from the disorder characteristic of 
uniform distributions. 

In addition to the above validations performed with 
synthetic data, the proposed approach for characteriza- 
tion of multiple-object systems was also applied to real 
data regarding the spatial distribution of the center of 
mass of the soma of retinal ganglion cells, kindly provided 
by H. Wassle. Reflecting the hexagonal organization of 
the photoreceptors characterizing the central portion of 
the mammals' retina, the position of these cells at the last 
retinal layer are also expected to follow some hexagonal 
organization, with some degree of disorder. The obtained 
rings around each original point and the respective func- 
tion h(s) are shown in Figures^a) and (b), respectively. 
The nearest-neighbor distribution considering uniformly 
spaced points for the same density are also shown (solid 



line) superposed to the curve h(s) in Figure It is clear 
from this result that the ganglion cell distribution is sub- 
stantially narrower than its respective uniform counter- 
part, corroborating a more organized scheme such as the 
hexagonal spacing. 

All in all, a new and effective approach to the 
characterization of multiple-object systems, based on 
artificially-induced percolation, has been proposed that 
presents a series of unique and interesting features. By 
considering only the exact distances representable in dis- 
crete spaces, as well as the fact that possibly co-existing 
point densities in such systems are typically sampled by 
using the same window, it has been possible to derive a 
powerful normalized tool capable not only of estimating 
the density of co-existing object distributions but also 
their respective perturbation degree. In addition, the 
number of points present in the observed region is re- 
flected by the height of the obtained peaks. It is argued 
that such a function provides one of the most informa- 
tive and accurate means for quantification of densities 
and spatial coverage, with the additional advantage that, 
unlike alternative schemes based on fractal and lacunar- 
ity, the obtained curve h(s) has an immediate and intu- 
itive interpretation in terms of point densities. Although 
not illustrated in this work, the proposed method can 
be immediately applied to multiple-object systems in- 
volving objects with shape and size more general than 
points. In such cases, the typical size of the objects has 
also been considered for the correct normalization and 
interpretation of the results. Other perspectives to be 
pursued further regard the consideration of the hierar- 
chical structures (a tree) obtained as a byproduct of the 
percolation-induced merging procedure and the possibil- 
ity of applying the artificially-induced percolation anal- 
ysis to characterize weighted networks. The proposed 
approach is currently being applied to characterize the 
texture of bones, growing neural systems, and spatial 
patterns of gene expression. 
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FIG. 2: Number of cluster mergings C(s) (a-d) and respective h(s) functions (e-h), represented by thick lines, obtained for 
uniform point densities with typical spatial scales 3, 7 11 and 15. The respective nearest-neighbor distributions are shown by 
the solid curves 
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FIG. 3: Profiles of h(s) (dotted lines) obtained for hexago- 
nal point distributions with increasing levels of perturbation. 
The nearest-neighbor distribution obtained for uniform point 
dispersion with equivalent density is shown by the solid curve. 
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FIG. 4: Spatial distribution of retinal ganglion cells (a) and 
respective h(s) function (b), represented by the dotted line. 
The nearest-neighbor distribution obtained for uniform point 
dispersion with equivalent density is shown by the solid line. 



